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q ! Abstract 
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Antivortices in ferromagnetic thin-film elements are in-plane magnetization configurations with 

Th" 1 ' 

' a core pointing perpendicular to the plane. By using micromagnetic simulations, we find that 

43 : 

magnetic antivortices gyrate on elliptical orbits similar to magnetic vortices when they are excited 
by alternating magnetic fields or by spin-polarized currents. The phase between high-frequency 

a ■ 

' excitation and antivortex gyration is investigated. In case of excitation by spin-polarized currents 

the phase is determined by the polarization of the antivortex, while for excitation by magnetic 

o : 

O . fields the phase depends on the polarization as well as on the in-plane magnetization. Simultaneous 

I \ excitation by a current and a magnetic field can lead to a maximum enhancement or to an entire 

\jq ' suppression of the amplitude of the core gyration, depending on the angle between excitation and 

m ! 

£NJ ■ in-plane magnetization. This variation of the amplitude can be used to experimentally distinguish 

\ between spin-torque and Oersted-field driven motion of an antivortex core. 
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I. INTRODUCTION 



Magnetic vortices and antivortices exist in ferromagnetic thin-film elements, where the 
interplay of demagnetization and exchange energy forces the magnetization out of plane to 
form a core in the centers.- 1 ^ The orientation of the vortex or antivortex core, denoted as the 
polarization p, is highly interesting for technical applications, e.g. magnetic memory devices, 
as it can be binary-coded.- 1 ^ Magnetic vortices have been studied intensively in the last years. 
It has been shown that a vortex core is deflected from its equilibrium position when excited 
by magnetic fields or spin-polarized currents.- 1 ^ The deflection causes a magnetic stray field 
which in turn exerts a force on the core.-^ The resulting gyroscopic motion can be described 
by a damped two-dimensional harmonic oscillator.- 

The dynamics of magnetic antivortices has hithertofore not been studied as intensively 
as magnetic vortex dynamics. Antivortices appear, e.g., in cross-tie domain walls and indi- 
vidual antivortices have been found in clover-shaped samples.-^'ii As illustrated in Fig. [H 
their in-plane magnetization shows a twofold rotational symmetry that is different from the 
continuous rotational symmetry of a vortex state. Due to their different in-plane magne- 
tizations, antivortex dynamics differs from vortex dynamics as is shown in this paper. An 
understanding of the dynamics of both, antivortices and vortices, is crucial for the descrip- 
tion of vortex-antivortex creation and annihilation. These processes have recently received 
a lot of attention as they are predominant features in the motion of cross-tie walls and in 
the switching of vortex cores.- 1 ^ 1 ^ 1 ^ 1 ^ 

Here we investigate the dynamics of antivortex cores, i.e. sense, phase, and amplitude of gy- 
ration, and compare them to the dynamics of magnetic vortices. We show that the direction 
of the in-plane magnetization around the (anti-)vortex core determines the phase between 
the exciting alternating magnetic field and the deflection of the (anti-)vortex core. For spin- 
polarized alternating currents the direction of the in-plane magnetization has no effect on 
the phase. Both micromagnetic simulations and an analytical model show that simultaneous 
excitation by magnetic fields and spin-polarized currents can lead to an enhancement or to 
an entire suppression of the antivortex core displacement. 

To classify vortices and antivortices the in-plane magnetization can be described by the 
relation^ 

= n(3 + 0o (1) 
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Figure 1: Definition of c = 2/ir ■ (<fi — n[3) for vortices and antivortices by Eq. (1). (a) Magnetic 
vortex (n = 1) with c = 1. (b) Antivortex (n = —1) with c = 1. (c) Antivortex (n = —1) with 
c = 2. 

between the angular coordinate of the local in-plane magnetization <fi and the angle (3 in 
real space with respect to the center of the (anti-)vortex core, as shown in Fig. [TJ The 
angles <p and (3 follow the mathematical sense of rotation. For a vortex n = 1 so that the 
in-plane magnetization turns in the same direction as the angle in real space with a constant 
difference 0o between and /3. For vortices, the angle 0o is independent of the choice of 
the axis to which (3 and <fi are measured. Thus for vortices <fro is an intrinsic quantity which 
can be expressed by the chirality c as 0o — ctt/2. In standard geometries and ferromagnetic 
materials stable vortices can only possess the chiralities c = 1 or c = -1. They can be 
mapped onto each other by mirroring the sample. In case of an antivortex n = — 1. This 
means that the in-plane magnetization turns opposite to the angle in real space. Though 
for antivortices the angle 0o is generally not conserved, because rotations of the sample lead 
to different values, we define a quantity c = 200/""" for antivortices with respect to a distinct 
axis.— Antivortices exhibit values c in the interval (—2,2]. A rotation of the antivortex 
by an anlge of leads to a change of the c-value of c = 20. This is due to the two-fold 
symmetry of the in-plane magnetization of an antivortex. 

II. MICROMAGNETIC SIMULATIONS 

To simulate magnetic-field induced antivortex dynamics the OOMMF— code sped up by 
higher order Runge-Kutta algorithms is used. The extended code includes the spin-torque 
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200 nm 500 nm 

Figure 2: Size and shape of (a) the vortex and (b) the antivortex sample. 

terms in the Landau-Lifshitz-Gilbert equation as given by Zhang and Li^-^, 

dM 



dt 



j'M x (H eS + x H eS 



- C 1 + a ^X M x ( M x U ■ V)M) (2) 

8 

-((-a)iMx(j.V)M, 

In this equation 7' = 7/(1 + a 2 ), where 7 is the gyromagnetic ratio, a the Gilbert damping, 
and £ the ratio between exchange and spin-flip relaxation time. The coupling between local 
current j and magnetization M is represented by = /xsP/[eM s (l + a 2 )], where P is the 
spin polarization. We simulate the excitation of a vortex in a 200 x 200 x 20 nm 3 permalloy 
square and an antivortex in a 500 x 500 x 40 nm 3 clover-shaped sample. The two geometries 
are shown in Fig. 2. Different thicknesses of t = 20 nm for the vortex and t = 40 nm 
for the antivortex sample are chosen in order to obtain similar eigenfrequencies for the 
two geometries. We assume a saturation magnetization M s = 8.6 • 10 5 A/m, an exchange 
constant A = 1.3 • 10" 11 J/m, a Gilbert damping parameter a = 0.01, and a ratio £ = 0.9a 
between exchange and spin-flip relaxation time.— ^ A lateral cell size of 4 nm is used. Thus 
the cell size is below the exchange length of permalloy of l ex = \J2A/ ^qM^ 5.3 nm. The 
position of the core is defined as the position of the maximum out-of-plane magnetization. To 
increase the spatial resolution, the magnetization of adjacent cells is matched by a polynomial 
of second order .— 
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Figure 3: Simulation of one gyration period of (a) an antivortex and (b) a vortex. Both have 
the topological charge q = —1/2 and are excited at the resonance frequency (727 MHz for the 
antivortex and 700 MHz for the vortex) by a current of amplitude j ■ P = 1.5 • 10 10 A/m 2 . The 
graphs below the magnetization images show the deflection in x-direction. (c) Exciting alternating 
current. 




III. EXCITATION BY MAGNETIC FIELD AND SPIN-POLARIZED CURRENT 



The eigenfrequencies of vortex and antivortex are determined by exciting the core with a 
current. The free relaxation of the magnetization yields the free frequency uj and damping 
T of the vortex or antivortex. Because the damping is small compared to the free frequency 
(r <C ujf), antivortex and vortex are weakly damped systems. Thus the frequency of the free 
oscillation uj and the resonance frequency u r are approximately the same. The simulated 
gyration of an antivortex core [p = 1, see Fig. 3(a)] and a vortex core [p = —1, see Fig. 3(b)] 
both driven by an ac current of amplitude j ■ P = 1.5 • 10 10 A/m 2 are shown in Fig. 3. 
They both possess the same sense of gyration as denned by the topological charge 2 ^ 2 ^ 2 ^ q = 
np/2 = —1/2. It is known that vortices gyrate counterclockwise with positive and clockwise 
with negative polarization.- Antivortices, on the other hand, gyrate clockwise with positive 
and counterclockwise with negative polarization. This is directly observed in the simulations. 
For small current or small magnetic field amplitudes, the simulated displacement of the 



antivortex is found to increase linearly with increasing excitation amplitude. This is due to 
the harmonic potential of the domains' stray field for small displacements of the antivortex 
core. Thus there is a linear restoring force on the antivortex, which has also been found for 
a vortex. 9 

To obtain the resonance curve for the amplitude and the phase t] of the antivortex core, 
either sinusoidal currents or magnetic fields of frequencies at, above, and below the resonance 
frequency u r are applied. Throughout this paper, the current is applied in x-direction while 
the magnetic field is applied in ^-direction. The resonance curve of a harmonic oscillator— 
with a resonance frequency uj r /2n = 727 MHz and a damping Y /2n = 6.4 MHz matches 
very well the numerical data, as shown in Fig. 4 (a,b). In general the antivortex gyrates 
on elliptical orbits. The semi-major (semi-minor) axis of the ellipses at frequencies below 
resonance changes into the semi-minor (semi-major) axis at frequencies above resonance. At 
resonance the trajectories are circular. This is illustrated in Fig. 4 (a) for a current-driven 
antivortex with c = 2 and in Fig. 4 (b) for a magnetic-field driven antivortex with c = 0. 
For both the semi-major axes point in y-direction at frequencies below and in x-direction at 
frequencies above resonance. 

The phase t] is defined by the temporal delay between the maximum of the applied current 
or field and the maximum core displacement in x-direction. Like for a harmonic oscillator 
the phase t] changes by ir when the exciting frequency Q is increased from values well below 
to values well above the resonance frequency u r . This is illustrated in Fig. 4 (a) for current 
and in Fig. 4 (b) for magnetic-field excitation. 

We numerically simulate the dependence of the phase on the direction of the in-plane 
magnetization by exciting at resonance antivortices of all possible integer c-values and of 
both polarizations p = — 1 and p — 1. For current excitation for all c-values the antivortex 
cores are deflected into the physical current direction (rj = ir). For magnetic field excitation, 
the phase is found to depend on the direction of the in-plane magnetization as shown in 
Fig. 5. For a constant frequency the phase varies by 2tt when c is changed from —2 to 2, i.e. 
when rotating the sample by n with respect to the magnetic field. At resonance the phase 
changes from to 2tt for p — 1 and from n to — n for p = — 1 as illustrated in Fig. 5. 

The phase r\ and its dependence on the direction of the in-plane magnetization is studied 
analytically by using the equation of motion for vortices and antivortices^ assuming low 
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Figure 4: (Color online) Exemplary resonance curves for the semiaxes of the elliptical trajectories 
and the phase r/ of an antivortex core gyration excited by (a) a current of amplitude j ■ P = 
1.5 x 10 10 A/m 2 (c = 2) and (b) a magnetic field of amplitude H = 300 A/m (c = 0). The 
symbols are results from micromagnetic simulations. The open squares illustrate the semiaxes in 
x-direction, the open circles the semiaxes in y-direction. The triangles show the phase r]. The 
asterisks illustrate the ratio between the semiaxes x and y. The solid line is the x-component and 
the dashed line is the y-component of the amplitude of a fitted resonance curve of a harmonic 
oscillator. The dotted red line is a fit of the phase i]. The insets show fits of the ratio between the 
semiaxes x and y as a function of the frequency. 
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Figure 5: (Color online) Phase rj between excitation and displacement for an antivortex core 
gyration at resonance. The dotted red line illustrates the phase when the antivortex core is excited 
by a spin-polarized current. The asterisks represent corresponding results from micromagnetic 
simulations for positive core polarization (p = 1) and the triangles for negative polarization (p = 
—1). For excitation with a magnetic field the solid line illustrates the phase for positive core 
polarization (p = 1), the dashed black line for negative polarization {p = — 1). The numerical 
results are depicted by open circles and squares. 

damping (uj r >> Y). The equation for the deflection 

= _ x . f v H sm(f)uj + (v H pcos(f ) + Vj ) iQ \ _ ^ iQt ^ 
V (v H ncos(^-) + Vjnp) uj — VffTip sin(^)zQ/ 

is derived from the Thiele equation^ for excitation with magnetic fields and the extension 
by Thiaville^ for spin-polarized currents. A harmonic potential due to the demagnetizing 
field is assumed.— The velocity due to the adiabatic spin-torque term is Vj = bjjo, the 
velocity due to the magnetic field vh = lHol/(27r), and the susceptibility of a harmonic 
oscillator is x — l/t^ 2 + + T) 2 ]- Equation (3) states that a change of the c-value leads 
to a rotation of the magnetic force. This in turn causes a dependence of the phase on the 
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in-plane magnetization. For example, a change from c = — 1 to c = 1 is equivalent to a 
rotation of the magnetic force by an angle of ir. At resonance Eqn. (3) yields the deflection 

x{c^t) = V ^e-^). (4) 

The maximum excitation is reached for iVtt = —m{ 1+p ~ pc ) ) so that the phase of the antivor- 
tex motion that is induced by a magnetic field at resonance reads 

fl+p-pc\ 

For purely current- driven excitation Eqn. (3) gives 

x(cp,t) = -%Le**. (6) 

Hence the phase induced by a current is rjj = it at resonance. This means that for current 
excitation the phase is independent of the direction of the in-plane magnetization as well as 
of the polarization. Figure demonstrates that the analytical results agree well with the 
simulations. 



IV. AMPLITUDE VARIATION OF GYRATION 



In the following we simulate antivortices that are excited simultaneously by a magnetic 
field and a spin-polarized current. First the antivortex is excited by a current in x-direction in 
the absence of a magnetic field. Then the amplitude of a magnetic field in ^-direction is tuned 
until the antivortex core gyration possesses the same amplitude as under current excitation. 
In our case the core is excited by a spin-polarized current of amplitude j ■ P = 0.7 TO 10 A/m 2 
that corresponds to a magnetic field of amplitude H = 122 A/m. Then the current and the 
magnetic field are applied simultaneously. Different directions of the in-plane magnetization 
(see Fig. 6) are chosen to investigate the c-dependent variation of the core amplitude for 
a positive polarization. The simulation shows a doubling of the amplitude at c = and 
a complete suppression at c = 2. Thus a superposition of the deflection by current and a 
perpendicular field leads to an amplitude variation in dependence on the direction of the in- 
plane magnetization of the sample. This is due to the c-dependent phase between antivortex 
core displacement and magnetic field, see Eqn. (131). The forces due to current and magnetic 
field are proportional to the deflections. If they are parallel or antiparallel, an enhancement 
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Figure 6: Amplitude of displacement of an antivortex core with polarization p = 1 at resonance. 
The antivortex is excited simultaneously by a magnetic field in y-direction and a current in x- 
direction. The symbols denote simulated results for vh = Vj, the lines are fits for different ratios 
vn/vj according to Eqn. (7) . 

or suppression of the core displacement is found, respectively. When both deflections have 
the same amplitude, the amplitde of gyration can be doubled or completely quenched, as 
shown in Fig. IIVI 

Using the addition theorem one can derive from Eqn. (3) the c-dependent amplitude variation 




of the antivortex core gyration at resonance. This is a general expression for arbitrary ratios 
vn/vj between the antivortex core velocities due to current and to magnetic field, for both 
polarizations and all c-values. The amplitude is plotted in Fig. 6. 

An inhomogeneous current distribution in the direction of the film normal generates 
a non-zero Oersted field perpendicular to the current.— For the experimental proof of a 
dependence of the amplitude on the direction of the in-plane magnetization, we propose 
a setup with a clover-shaped sample that is illustrated in Fig. [7J A similar sample was 
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Figure 7: Proposed setup with electrical contacts to excite a single antivortex, here with polariza- 
tion p = 1. The quantity c depends on the direction of the exciting ac current. For current in x- 
or y-direction the deflection is suppressed or amplified, respectively. 

investigated by Shigeto et al.— with magnetic-force microscopy. Excitation by a spin current 
in x-direction or y-direction through electrical contacts corresponds to a direction of the in- 
plane magnetization for c = 2 and c = 0, respectively. For the polarization p = 1 we expect 
a suppressed motion when the current is applied in x-direction and an enhanced amplitude 
when it is applied in y-direction. The variation of the amplitude for c = and c = 2 could 
be used to determine the ratio between the forces on the antivortex core due to an Oersted 
field and a current. 



V. CONCLUSION 



In conclusion, we have demonstrated by micromagnetic simulations that antivortices ex- 
cited by spin-polarized ac currents or magnetic fields gyrate on elliptical orbits. These orbits 
can be well described by the analytical model of a two-dimensional harmonic oscillator. The 
sense of gyration of antivortices depends solely on the topological charge q = np/2. The 
phase of the antivortex motion excited by an alternating magnetic field depends also on the 
direction of the in-plane magnetization. Antivortices that are excited simultaneously by a 
spin-polarized current and a magnetic field show an enhancement or a suppression of the 
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deflections amplitude in dependence on the direction of the in-plane magnetization. The 
effect of the amplitude variation in dependence on the in-plane magnetization can be used 
to experimentally investigate the influence of Oersted fields in current-induced antivortex 
dynamics. 
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